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The critical properties of the frustrated Kondo necklace model with a half saturation magneti- 
zation (m = 1/2) have been studied by means of an exact-diagonalization method. It is shown 
from bosonization technique that the model can be effectively expressed as a quantum sine-Gordom 
model. Thus it may show three (dimer plateau, Neel plateau and Tomonaga-Luttinger liquid) phases 
due to competitions among the Ising anisotropy A, and the nearest- and next-nearest-neighbor ex- 
change interactions Ji and J2. The boundary lines on the A — J2/J1 phase diagram separating 
the three phases are determined by the method of level spectroscopy based on the conformal field 
theory. 
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I. INTRODUCTION 

Intermetallic compounds containing rare-earth or ac- 
tinide elements are known to exhibit a variety of low- 
temperature phases, including antifcrromagnetism, su- 
perconductivity, and Kondo-insulator phases. The 
Kondo lattice model is one of the standard models de- 
scribing such physics of heavy electron systems. The 
low-temperature properties of the Kondo lattice model 
at half-filling, describing a group of compounds called 
Kondo- insulators, are governed by competition between 
the Kondo- and the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactions. While the Kondo interaction fa- 
vors the formation of intra-atomic quantum disordered 
(Kondo singlet) phase, the RKKY interaction favors an 
antiferromagnetic long-range order. In one-dimension, 
the Kondo singlet phase is stabilized in the presence of 
Kondo interaction of any strength, causing an energy gap 
in the charge and spin excitations, the spin gap being al- 
ways smaller than the charge gap.0, ^ 

The Kondo necklace model is a simplified version of 
the Kondo lattice model in which the charge degree of 
freedom is frozen out in order to describe the above- 
mentioned competition. || In fact, the Kondo singlet gap 
Ak always opens in the Kondo necklace model as long 
as any finite Kondo interaction is present. ||, |5|, |[ [?], || 
Thus the magnetization of this model starts to grow at 
a finite critical magnetic field h = Ak , and continuously 
increases with increase in the external field. Note that 
these results for the Kondo necklace model are consis- 
tent with the magnetization curve of the Kondo lattice 
model. H 

However, the above situation may be drastically 
changed when the next-nearest-neighboring conduction 
electron spins coupled antiferromagnetically. In our pre- 
vious paper, pp| we studied a magnetization process of 
the Kondo necklace model with next-nearest-neighbor 
interaction J2 which partially frustrates the nearest- 
neighbor interaction J± , making use of a finite size scaling 
analysis of the numerical exact-diagonalization data, 

In Ref. 10 we showed that for A = J2/J1 = 1/2, the 



magnetization plateaus appear at zero (m = 0) and a 
half saturation magnetization (to = 1/2), satisfying the 
necessary condition for the appearance of plateaus, Dj] 
The origin of these plateaus for A = 1/2 can be under- 
stood as follows. The spin configuration in the plateau 
phase at to = is interpreted as the aggregate of intra- 
atomic Kondo singlets. On the other hand, the spin con- 
figuration in the plateau phase at to = 1/2 represents 
the aggregate of intrachain dimers with broken transla- 
tional symmetry, resulting from the frustrations due to 
competition between J\ and J2. However, the to = 1/2 
plateau may be completely broken by quantum fluctua- 
tions with decreasing J2 In fact, the to = 1/2 plateau 
does not appear at a magnetization curve of the Kondo 
necklace model without next-nearest-neighbor interac- 
tion. This suggests that a quantum phase transition be- 
tween the plateau and nonplateau phases may occur in 
0<A<l/2. 

In this paper, we are going to investigate critical prop- 
erties of frustrated Kondo necklace model 

L 

^FKN = ^ {(Oj • Oj+i)a +M°J " °J+2)a) 
J=l 

L 

+J K Y,"i-Si, (!) 
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removing the restriction A = 1/2 imposed in our previous 
work.JlO| Here (a-j -crj + s) A with S = 1,2 denotes the 
anisotropic spin interaction of the XXZ type, 



((Tj ■ cr j+s ) A = a'a 



j+S 



Aa]a] +& 



(2) 



with Ising anisotropy A, crj are intra-chain spins with 
(7 = 1/2, and Sj are innercore spins with S = 1/2, 
hanging from the a-spin chain. We choose the nearest- 
neighbor exchange J\ as the unit of energy, and the 
Kondo exchange coupling Jk is chosen to be antiferro- 
magnetic (Jk > 0). In the case of the ferromagnetic 
Kondo coupling (Jk < 0), this model ([!]) has a connec- 
tion with the Haldane problem.) 12 
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We are also interested in whether or not the above- 
mentioned non-plateau phase belongs to the universality 
class of Tomonaga-Luttinger (TL) liquid. This problem 
is nontrivial in the case of the Kondo necklace model, al- 
though it is well known that many one-dimensional quan- 
tum spin systems exhibit the TL liquid phase with central 
charge c = 1 in the conformal field theory. 

The paper is organized as follows. In we map 
model ([TJ) onto the quantum sine-Gordon model by 
means of bosonization technique, limiting ourselves to 
the Hilbcrt space with the half-magnetized states. In 



III A, the phase diagram in the plane of Ising anisotropy 



A and the frustration parameter A will be constructed, 
making use of the method of level spectroscopy, on the 
basis of the assumption that our model (IH) is described as 



the quantum sine- Gordon model. In STUB, the validity 



of this assumption will be checked by calculating scaling 
dimensions and the central charge. It will be shown that 
the central charge c is very close to 1. The final section 
(§ IV) is devoted to a summary and discussion. 



II. BOSONIZATION AND EFFECTIVE 
HAMILTONIAN 

Let us consider the effective Hamiltonian that de- 
scribes our model in the m = 1/2 plateau phase. We 
showed in our previous paper ]l0| that inner core spins of 
the Kondo necklace model with A = 1/2 in the to = 1/2 
plateau phase are almost ferromagnetically aligned along 
the external field (z-axis) . This enables us to replace the 
inner core spin operators Sj with an average (Sj) = 1/2, 
leading to the effective Hamiltonian of the following form: 

L 

Ti-cs = 22{(crj ■ cr j+ i) A + A (a j ■ <T j+2 ) A } 

3=1 



3=1 



(3) 



Here h c g — (h — Jk/2) is the effective magnetic field in- 
cluding the effect of Kondo coupling. Since the total mag- 
netization of the a spins is vanishing (Jr-/2 (h e s ~ 0), 
the system can be effectively regarded as a a — 1/2 XXZ 
chain with next-nearest-neighbor interaction in the zero 
field. The ground state properties of the model have been 
extensively studied by many authors by means of various 
methods . |L3|, [l4|, [l5| [l(| It is known that the model is 
transformed into the quantum sine-Gordon model 



dx 



(d x (j)(x)) 2 



2-kK 



2.9 



(27ra) z 
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in terms of the Jordan- Wigner transformation and 
bosonization technique. Here, n(x) is the momentum 
density conjugate to the bosonic field 4>{x), satisfying the 



commutation relation [<j)(x), IL(x')] — i8(x— x'), and a is 
the lattice constant. The spin wave velocity v s , the Lut- 
tinger liquid parameter K, and the umklapp scattering 
amplitude <? are respectively given as 



K = 



g = 2tt 2 \ 2 D, 



(5) 



where the coefficients A, C, and D are given by 
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Thus the critical properties of our system with to = 1/2 
are well described in terms of the quantum sine-Gordon 
model, although the above expressions for A, C, and D 
are valid only when A, A and h e ff are small. It is well 
known that the quantum sine-Gordon model exhibits two 
gapful (corresponding to dimer plateau and Neel plateau) 
phases and a gapless (TL liquid) phases, depending on 
the value of K. These results are derived from the flow 
diagram of renormalization group equations derived from 
the scaling of the cut-off a — > e d 'a where I = InL. 

Thus the frustrated Kondo necklace model belongs to 
a universality class of TL liquid if the nonlinear term is 
renormalizcd to zero as L — > oo. In the TL phase the 
CT-spin excitation is gapless, and the cr-spin correlation 
functions algebraically decay as 



A" 



(-iy> X)-<O^ 



(7) 

(8) 



On the other hand, the dimer plateau phase is charac- 
terized by the cr-spin excitation gap, the exponential de- 
cay of the c-spin correlation, and dimer long-rang or- 
der. The Neel plateau phase is characterized by the 
Ising gap for cr-spin and antiferromagnetic long-rang or- 
der (Neel order). Also, it is known that the dimer-TL and 
the Neel-TL transitions are of the Berezinskii-Kosterlitz- 
Thouless (BKT) type, and the dimer-Neel transition is 
of the Gaussian type. 



III. NUMERICAL RESULTS 

A. Phase diagram 

In this section, we are going to determine the bound- 
aries between TL, dimer, and Neel phases on the plane 
of Ising anisotropy A and the frustration parameter A, 
making use of the method of level spectroscopy. This 
method is known to be powerful in the determination 
of the BKT transition points, which is difficult by the 
standard finite size scaling analysis due to logarithmic 
corrections. |0| [H], [TJ] 
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TABLE I: Symmetry classification of the ground state and 
three (spin-wave type, dimer and Neel excitations) excita- 
tions in the Kondo necklace model with m — 1/2 and L — 
4 x (integer). The numbers in round brackets are quantum 



numbers in case of L = 


4 x integ 


cr + 2. 










Quantum Numbers 








M 


k 


V 




ground state 


L/2 


(tt) 


+1 (- 




spin-wave excited state 


L/2± 1 


7T (0) 


-1 (H 


-1 


dimer excited state 


L/2 


7T (0) 


+1 (" 




Neel excited state 


L/2 


7T (0) 


-1 (H 





Since the Hamiltonian (|l|) is invariant under spin rota- 
tion around the z-axis, translation j — > j + 1, and space 
inversion j — ► L — j + 1, the eigenvalues of total magneti- 
zation M, wave number k = 2im/L, and parity V = ±1 
are good quantum numbers to label the eigenvalues and 
eigenvectors of Hfkn- 

As expected from bosonization, the ground state with 
to = M/L = 1/2 in all three phases (i.e., TL, dimer 
and Neel phases) is shown to be labeled by (M, k,V) — 
(L/2, 0,1) for the lattice with L = 4 x (integer) and 
(M,k,V) = (L/2,71-,-1) for the lattice with L = 4 x 
(integer) + 2 On the other hand, the lowest excited state 
in TL, dimer, and Neel phases are labeled by (M, k, V) = 
(L/2 ± 1,tt,-1), (L/2,t:, 1), and (L/2,tt,-1) for L = 
4 x (integer) and (M,k,V) = (£/2± 1,0,1), (L/2, 0,-1), 
and (L/2, 0, 1) for the lattice with L = 4 x (integer) + 2, 
respectively. Hereafter, we call each of the excitations in 
TL, dimer, and Neel phases spin-wave type, dimer, and 
Neel excitaions. Table | summaries the eigenvalues of 
these symmetry operators in the ground state and three 
(spin-wave type, dimer and Neel excitations) excitations 
in the frustrated Kondo necklace model with to = 1/2. 
In the case of L = 4 x (integer), the (finite size) critical 
points are determined from the crossing of the minimum 
excitation energies of the above three types of excitations 

AE SW = \{E(M = L/2+l,k = TT,P= -1) 

+E(M = L/2 - 1, k = tt,P = -1) 

-2E(M = L/2,k = 0,P = 1)|, (9) 

A£ dimcr = E(M = L/2, k = TT,P=l) 

-E(M = L/2,k = 0,P = 1), (10) 

A^ N6 el = E(M = L/2,k = TT,P=-l) 

-E(M = L/2,k = 0,P = 1). (11) 

Note that these minimum excitation energies corre- 
spond to the scaling dimensions x sw , Xdimcr, and XNcci 
by 

Xi = — AEi (i — sw, dimer and Neel), (12) 



where v s is the spin wave velocity defined as 

v s = lim — { E (M = L/2, k = 2n/L) 

-E(M = L/2,k = Q)\. (13) 

We defined the gap in the spin-wave type excitation 
AE SW as in eq.(||) so that the Zeeman energy due to the 
external magnetic field 

L 

H Zccman = -hY,{°j+S*) (14) 

i=i 

is canceled out. 

Let us consider the TL-dimer critical point for A = 
Jfc = 1.0 in order to illustrate how our method of level- 
crossing works. In Fig.|l|, spin-wave, dimer and Neel ex- 
citation gaps for L = 8 are shown as functions of the 
parameter A. We see from Fig.|l| that the spin-wave type 
excitation is the lowest for A < 0.35105 and so is the 
dimer excitation for A > 0.35105. On the other hand, 
the Neel excitation is always higher than spin-wave type 
and/or dimer excitations. 

We have estimated the values of A at the crossing 
point of spin-wave type and dimer excitations for vari- 
ous sizes L = 6,8, 10, 12 and 14. These are extrapolated 
to the thermodynamic limit L — > oo, assuming a formula 
X C (L) = \ c + A/L 2 + B/L 4 . Thus the TL-dimer critical 
point A c = 0.3419 ±0.0001 (See Fig.|) is obtained in the 
limit L — > oo. We have performed similar procedures 
for various A's to construct the phase diagram on A — A 
plane, as shown in Fig|| We assumed that to = 1/2 and 
Jk = J\ in these calculations. 



° dimer 




FIG. 1: The dimer, Neel, and spin- wave type excitation ener- 
gies as functions of the frustration parameter A. The crossing 
point of dimer and spin-wave type excitaion energies repre- 
sents the BKT transition point for L = 8 and A = 1.0. 
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FIG. 2: The crossing points of dimer and spin- wave type exci- 
tation energies for various sizes (L = 6, 8, 10, 12 and 14) and 
A = 1.0. In the thermodynamic limit, the BKT transition 
point is obtained as A c = 0.3419 ± 0.0001, assuming an ex- 
trapolation formula A C (L) = A c + A/L 2 + B/L 4 , assuming an 
extrapolation formula X C (L) = A c + A/L 2 + B/L 4 . 
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FIG. 3: Phase diagram on the A — A plane of the frustrated 
Kondo necklace model with m — 1/2 for ,1k = 1.0. The open 
circle denotes a multicritical point between the dimer, Neel 
and TL phases. 



B. Consistency Check 

In §|0|, we assumed that the effective Hamiltonian of 
our model at m = 1/2 is described by the quantum sine- 
Gordon model. Our analysis, making use of the method 
of level spectroscopy, is based on this assumption, accord- 
ing to which the scaling dimension without logarithmic 
corrections, 



FIG. 4: The average scaling dimension x — (xdima + ^Neei + 
2:e S w)/2 along the TL-dimer and dimer-Neel boundaries in 
Figg 



TABLE II: The relations between the scaling dimension 
[Eqs. © and (0)]. 



A 


^Nccl/ 2-dimcr 


^sw(^dimcr H~~ ^Nccl)/2 


2.6 


0.999 


0.2492 


2.8 


1.004 


0.2490 


3.0 


1.005 


0.2486 



should tend to 1/2 on the BKT critical line in the ther- 
modynamic limit L — > oo. We confirmed this relation 
x w 1/2 on the TL-dimer abd dimer-Neel transition lines, 
as clearly shown in Fig.|]. Furthermore, Table. [n] indicates 
that the relations between the scaling dimensions 



^Necl — ^dimcr, 
^sw-^dimcr 1/^ 



(16) 
(17) 



are satisfied on the dimer-Neel (Gaussian) transition 
points. 

We have also estimated the central charge c along the 
TL-dimer and dimer-Neel critical lines from the finite size 
correction to the ground-state energy 



6L 



(18) 



where eo is the ground-state energy per site (including 
two spins) in the thermodynamic limit. As shown in 
Fig.||, we obtained the value c w 1 with small errors. 
These results provide sufficient evidence that our model 
can be described by the quantum sine-Gordon model. 



IV. CONCLUSION 



^dimer *^Neel 2x s 



In this paper, we have studied critical properties of the 

(15) 

frustrated Kondo necklace model with a half-saturation 
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ing dimension without logarithmic corrections [Eq. JTe|)] 
are close to 1/2 at the TL-dimer (BKT) and dimer-Neel 
(Gaussian) transition points. The relations between the 
scaling dimensions [Eqs. ([l6]) and (|l7|)1 are satisfied as 
the Gaussian transition points. Furthermore, the central 
charge c ~ 1 is obtained along the BKT and Gaussian 
transition lines, by estimating the finite-size correction 
for the ground state energy [Eq. (pq)]. Thus we can 
conclude that the present system belongs to the same 
universality class as the quantum sine-Gordon model. 
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FIG. 5: The central charge c on the TL-dimer and Neel-dimer 
critical lines, yielding the result c = 1 within errors of 0.03%. 



magnetization (m = 1/2), making use of the method of 
level spectroscopy. The phase diagram on the plane of 
the Ising anisotropy A and the frustration parameter A is 
constructed exhibiting two plateau (dimer and Neel) and 
a TL phases. We also confirmed that the values of scal- 
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